Attosecond spectroscopy of solids: streaking phase shift due to lattice scattering 
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Theory of laser-assisted photoemission from solids is developed for a numerically exactly solvable 
model with full inclusion of band structure effects. The strong lattice scattering in the vicinity of 
band gaps leads to a distortion and a temporal shift of the streaking spectrogram of the order of 
100 as. The effect is explained in terms of Bloch electron dynamics and is shown to remain large for 
an arbitrarily small photoelectron mean free path. The implications for the streaking experiment 
on W(110) [A. L. Cavalieri et at, Nature (London) 449, 1029 (2007)] are discussed. 

PACS numbers: 79.60.-i, 34.80.Qb, 42.65. Re 



Attosecond photoelectron streaking is a unique tool to 
study ultrafast electronic processes with subfemtosecond 
resolution [l[ . The method relies on mapping time to en- 
ergy by means of a strong laser field: the electron wave 
packet created by an ultrashort pulse of extreme ultravi- 
olet radiation (XUV) is accelerated by the superimposed 
laser field, and its energy spectrum J(e) shifts up or down 
depending on the momentum transfer from the laser field 
to the outgoing electron. The energy gain depends on the 
electron release time r relative to the laser field E h (t), as 
encoded in the spectrogram J(e, r). In comparison to the 
isolated atom, the application of the method to solids is 
complicated by two factors. First, the outgoing wave 
is a coherent superposition of contributions from many 
atoms. Second, the crystal lattice strongly affects the 
propagation of the photoelectron. 

In the seminal experiment on W(110) by Cavalieri et 
al. a shift of Ar = 110 ± 70 as was discovered be- 
tween the spectrograms of valence 5d and semi-core 4/ 
states. Its origin has been debated in recent theoretical 
studies In Refs. 0-0] the original idea Q is devel- 

oped that the shift is due to a spatial inhomogeneity of 
the streaking field: under the assumption that the laser 
field does not penetrate into the crystal the measured At 
is just the difference in the travel time to the surface for 
bd and 4/ electrons. Baggesen and Madsen proposed 
to use the states located outside the metal as a reference 
to measure the absolute travel times. In the classical 
simulation by Lemell et al. 4| the streaking phase shift 
was traced back to the escape depth and to the electron 
group velocity. Kazansky and Echenique [5], however, 
questioned the applicability of the concept of group ve- 
locity as derived from band structure and ascribed the 
delay both to the different character of the initial 5d 
and 4/ states and to a different propagation of the out- 
going electrons. Contrary to the other authors, Zhang 
and Thumm [6j assumed the laser field to be spatially 
constant up to the escape depth of 5 A, but, similar to 
Ref. @, found the difference in the localization of the 
initial states to be the reason for their different streaking 
behavior. One important aspect not considered in previ- 



ous studies is the Bloch nature of the outgoing electron 
and its implications for the streaking spectrogram. 

Knowledge of the crystal lattice scattering effects is es- 
sential for the understanding of what is being measured in 
a streaking experiment on a solid. The aim of the present 
Letter is to shed light on how the band structure affects 
the formation of the spectrogram. Its role is twofold: 
First, close to the band gaps, the group velocity of the 
outgoing electron is reduced, which reduces the energy 
transfer to the photoelectron even if the laser field in the 
solid is as strong as in vacuum. Second, if the final state 
is close to a band gap edge, an electron at the surface 
is brought by the XUV pump pulse towards the crystal 
interior before it goes back and leaves the solid. This 
complicated trajectory causes a distortion of the spec- 
trogram, which can be perceived as a final-state-energy 
dependent delay in photoemission. The joint effect of the 
two factors is shown to result in values of Ar of the order 
of 100 as, i.e., close to those experimentally observed. 

To arrive at reliable conclusions on the dependence of 
Ar on the final state energy it is necessary to avoid any 
uncontrollable numerical approximations. We consider 
a model crystal surface without lateral corrugation, for 
which the time-dependent Schrodinger equation (TDSE) 
can be solved and the spectrogram J(e, r) can be ob- 
tained with any desired accuracy. The periodic lattice is 
represented by a slab of 24 atomic layers with a piece- 
wise constant potential V(z), see Fig. QJa). We consider 
normal emission and assume a linear polarization of the 
light with electric field along the surface normal, so the 
problem reduces to a one-dimensional (ID) TDSE. The 
perturbing field E(t) is a superposition of an XUV and a 
laser pulse, see Fig.[TJc). The temporal envelopes of both 
pulses are of the form cos 2 '(irt/D), The full duration of 
the XUV pulse is D = 500 as. The laser pulse is chosen to 
be an odd function of t: E L (t) = E^ sin(fii) cos 2 (nt / 'D) , 
with photon energy 57 = 1.65 eV, D — 5 fs, and ampli- 
tude E™ = 2 x 10 7 V/cm. The syste m is enclosed in a box 
[Fig. Hlb)], an d all calculations are performed in matrix 
form in terms of exact eigenfunctions (discrete spectrum) 
of the unperturbed Hamiltonian H = p 2 /2m + V(z). The 
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FIG. 1: (color online) (a) Potential profile of the model crys- 
tal. Atomic layers are given by potential wells, lattice con- 
stant is 4 a.u. (b) Domain of H: —750 < z < 750 a.u. (c) 
Sum of XUV and laser pulses, (d) z-resolved density of states 
(DOS). Energy is relative to the "muffin-tin" (MT) zero. 



maximal at r = 0, but because at that point the elec- 
tric field increases the electron in the crystal gains more 
energy if it starts later to be exposed to a stronger field 
while it moves faster. Figure Ufa) shows the exact result 
for a classical particle (full line) starting from the depth 
of 55 a.u. with the initial kinetic energy of 94 eV in a 
solid with an inner potential of U = —30 eV. This corre- 
sponds to the excitation of the narrow band by the XUV 
pulse of w = 106 eV (e = 69 eV), and the fully quantum- 
mechanical calculation yields an almost identical streak- 
ing curve: it is seen to be shifted to the right by ~ 50 as 
from the free-particle curve @ . What happens when the 
final state is close to a band gap, where the electron group 
velocity in the crystal is lower than in vacuum? Here, by 
the Bloch electron dynamics law hdk/dt = F, its energy 
changes slower than in vacuum, and the opposite effect is 
expected. Indeed, for u = 155 eV (e = 118 eV) the s(r) 
curve is shifted to the left, see Fig.^b), i.e., the electron 
appears to be released earlier than the free particle. 



electric field perturbation is taken in the dipole approxi- 
mation, AH = zeE(t). Note that the XUV pulse is not 
separated from the laser pulse - input into AH is their 
sum. The TDSE is solved with the split-operator tech- 
nique in the spectral representation, so the effect of the 
crystal potential is fully taken into account for both ini- 
tial and final states. The spectrum of H is truncated at 
200 eV above the vacuum level, which ensures the conver- 
gence of all the results. More details of the methodology 
can be found in Ref. [7| , where the model was applied to 
a single atom. 

Initial states are the two bound bands, see Fig. [TJd): 
the narrow band (1) and the wide band (2), respectively, 
model semi-core and valence states. Each initial state 
gives rise to two wave packets traveling in opposite direc- 
tions. When both pulses are over and the right-traveling 
packet has left the crystal, it is reexpanded in terms of 
the eigenfunctions of H to yield the recorded spectrum. 
The observed spectrum J(e, r) is a sum of the spectra of 
the 24 levels representing the band. The centroid (center 
of gravity) of the total streaked spectrum e(r) shifts rela- 
tive to the centroid of the laser- free spectrum en, and the 
dependence of the shift on the release time r is referred 
to as the streaking curve, s(r) = e(r) — en. Note that 
for a classical particle created at the time r and moving 
under the field E L (t) in a flat potential V(z) = const, the 
streaking curve is symmetric: s(r) = s(— r). 

Let us start with the obvious observation that the aver- 
age potential in the crystal U is lower than the potential 
in vacuum, so in the crystal a classical electron moves 
faster than in vacuum. Thereby its streaking curve devi- 
ates from that of a free particle, and this distortion can 
be perceived as a delay of the release of the photoelectron 
coming from the depth of the crystal relative to the one 
released just at the surface. For example, in the present 
case, the momentum transfer for a free electron would be 
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FIG. 2: (color online) (a) Blue solid line is for a classical par- 
ticle traversing the potential barrier U . (a) and (b) Circles 
are the s(t) values by TDSE for the emission from band 1 
(e ini = -37 eV) with u = 106 eV (a) and 155 eV (b). Red 
dashed lines are for a classical particle (starting with the same 
kinetic energy) in a flat potential. Inset illustrates the defi- 
nition of a. (c) and (d) Streaking curve shift <r(e) for a clas- 
sical particle (thick dashed line) and for emission from band 
1 (black full circles) and band 2 (etni = —9 ± 5 eV, red open 
circles), e = eini+w. Laser field is E" = 2 x 10 7 V/cm (c) and 
4 x 10 7 V/cm (d). Shown to the right are the crystal band 
structure e(k) and the parabolic fit to the unbound spectrum 
(red dashed line) to determine the inner potential U. 
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For the special case of an odd function E h (t) the tem- 
poral shift of the curve can be roughly quantified as the 
center of gravity a of the positive part of the function 
s(t). Its dependence on the final energy er(e) for the clas- 
sical particle and for both bands is shown in Figs, ^c) 
and [Hd) ■ Far from the band gaps the streaking curve 
shows a delay, and a is close to that for a classical par- 
ticle starting at —55 a.u., which is slightly deeper than 
the average spatial location of the initial states, —48 a.u. 
When the final state energy approaches a band gap the 
s(r) curve strongly deviates from the free-electron shape, 
and the function er(e) shows oscillations. This is caused 
by rapid energy variations of the photoemission cross- 
section in combination with the large energy width of 
the pump pulse [9| , so that the spectrum is composed of 
several contributions, which are differently streaked by 
the probe field. Nevertheless, when the spectral weight 
concentrates close to a band-gap edge, the s(t) curve re- 
tains an approximately symmetric shape and shifts to the 
left, showing an advance (negative a). The er(e) curves 
are rather stable to the strength of the streaking field: a 
two times larger amplitude E™ yields almost the same re- 
sult [cf. Figs. [He) andj^d)]. This proves the robustness 
of the effect and that the parameter a is a reasonable 
choice to characterize its range. The <r(e) curves are very 
similar for the emission from the narrow and from the 
wide band, with more oscillations in the latter case due 
to the much broader spectrum of the wide band. Far 
from the gaps the values of a for the two bands closely 
agree, demonstrating the dominant role of final states. 

We may conclude that the semi-classical picture based 
on the notions of inner potential and group velocity is 
a good starting point to understand the r-shift of the 
streaking curve. To relate it to the experiment, we must 
take into account the surface sensitivity of photoemission: 
owing to inelastic scattering, only the electrons excited in 
a close vicinity of the surface contribute to the observed 
spectrum. The travel path to leave the crystal is, thus, 
limited by the escape depth A, and for sufficiently small 
A the effect of delay that occurs in nearly-free-electron 
(NFE) regions (far from band gaps) must become negligi- 
ble. However, the effect of advance would not completely 
disappear even for a small A of few Angstroms because 
close to the gap the velocity may be arbitrarily small. 
For an unambiguous analysis of the escape depth effect 
it is desirable to avoid a phenomenological treatment of 
inelastic effects by an imaginary absorbing potential and 
current non-conservation. To introduce a spatial resolu- 
tion into the microscopic theory, let us create a defect in 
the ID lattice and study the streaking of the photoelec- 
tron from the localized state at the defect. 

Figure [3] shows the energy dependent distortion of the 
streaking curves for the initial state localized at the defect 
in the first, second, and fourth layer. In accord with the 
above discussion, the effect of delay in the NFE regions is 
rather small (although it increases with the depth) . How- 
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FIG. 3: (color online) (a) and (b) Delay <r(e) for emission from 
the state at the defect (inset) in the 1st (black full circles), 2nd 
(red open circles), and 4th layer (blue crosses), eini = —41 eV, 
e = eini + uj- (c) and (d) s(r) for uj — 82 eV for the defect in 
the 1st layer (full lines) and for a classical particle in a flat 
potential (dashed lines). Vertical lines show band gaps, (a) 
and (c) are for the laser pulse E h (t)\ (b) and (d) for —E L (t). 



ever, the effect of advance is quite large: as in the case of 
the extended states, the variations of a are comparable 
to the measured At value. The most surprising is that 
for the 1st layer a does not vanish - it is of the same 
order or even larger than for the 2nd and the 4th layer. 
Although for the electron released at the 1st layer there 
seems to be no space to traverse to leave the crystal, it 
appears to have spent some time in the crystal. In fact, 
the wave packet is brought inside the crystal by the pump 
pulse. This happens when the final state is at a band-gap 
edge, where the local DOS is larger in the crystal than in 
vacuum, see Fig. Hfd). Owing to the DOS gradient, the 
wave packet is shifted beneath the surface. Comparison 
of the streaking by the pulse in FigQJc) and by the pulse 
of the opposite sign (acceleration at r = vs. decelera- 
tion) again indicates that the effect is a property of the 
system rather than of the specific perturbation. 

Some features of the <r(e) curves can be readily under- 
stood within the classical picture: the minima B and E 
and the maximum C are more pronounced for the deeper 
defects. The minima A and D, however, show the oppo- 
site trend. This happens when a rapid energy dependence 
of the photoemission cross-section leads to oscillations of 
a. Close to the surface this energy dependence varies 
with the position of the defect, and at certain energies 
the effect becomes stronger for the 1st layer than for the 
deeper layers. Nevertheless, also in this case the streak- 
ing curves retain their traditional shape; an example for 
the topmost layer is shown in Figs.[3][c) and[3fd). 

These results show that the streaking phase shift may 
occur at arbitrarily small escape depths and without the 
interference between the layers (in the case of the defect 
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FIG. 4: (color online) (a) Real band structure (thin lines) 
and conducting complex band structure (red thick lines) of 
W(110). (b) Energy dependence of group velocity for four 
conducting branches (full lines) and in vacuum (dashed line). 

it is explicitly excluded). This imposes a limitation on 
using the states located outside the metal as a zero-delay 
reference: the final states must be in an NFE region. 

Finally, let us see whether the present theory is rele- 
vant to the experiment on tungsten. The band structure 
for the normal emission from the (110) surface (NTN 
line) is shown in Fig. 2) To determine the Bloch waves 
that effect the electron escape into vacuum we draw on 
the one-step theory of photocmission, in which the fi- 
nal states are time reversed LEED states (low energy 
electron diffraction). The LEED wave functions are cal- 
culated from the complex band structure of the semi- 
infinite W(110) crystal with the linear augmented plane 
wave method as explained in Ref. . The Bloch waves 
most strongly contributing to the LEED state - the con- 
ducting branches - are highlighted by the line thickness, 
which is proportional to the current carried by the in- 
dividual partial wave (see Ref. [ll[ for the explanation). 
The energy dependence of the group velocity for the four 
conducting branches is shown in Fig. S{b). The final 
states for photocmission from both 4/ and 5c? bands are 
seen to fall in the vicinity of the band gaps, at 58 and 
91 eV, respectively. Thus, the effect of advance must 
occur in both cases, and, judging by the variations of a 
given by the present model (as large as 150 as), it strongly 
contributes to the observed effect. Qualitatively, one may 
speculate that because the velocity mismatch is larger at 
higher energies the negative shift of the streaking curve 
should be larger for 5d than for 4/ emission, which would 
give an advance of the 5d emission, as in the experiment. 

To summarize, the final states band structure causes 
a r-shift of the spectrogram, and its sign depends on 
whether the photoelectron moves in the crystal faster or 
slower than in vacuum. For an electron coming from a 
large depth and moving faster than in vacuum the streak- 
ing curve shifts towards positive delays. In the present 
ID model (consistent with realistic k|| projected band 
structure) the delay is within 50 as. This effect is in ac- 



cordance with classical dynamics. Naturally, it steadily 
decreases with decreasing escape depth. The quantum 
nature of the outgoing electron becomes important for 
final energies sufficiently close to a band gap, where it 
always moves slower than in vacuum. In that case the 
spectrogram shows an advance. The quantum effect is 
found to exceed 100 as even for electrons excited from the 
states localized at the surface. This reveals an aspect im- 
portant for various electron spectroscopies: the outgoing 
electrons are sensitive to the substrate band structure no 
matter how small the electron mean free path is. 

The main implication for the attosecond spectroscopy 
is that neither the inhomogeneity of the streaking field 
nor an extended initial state are a prerequisite for an ap- 
preciable relative delay. Still, the mean free path plays 
a large role because depending on its value different as- 
pects of the mechanism become important. Experimen- 
tally, the predicted effect should manifest itself as a non- 
monotonic dependence of At on the XUV energy, and it 
must be allowed for in interpreting spectrograms. At the 
same time, measurements at different uj would be instru- 
mental in separating the final-states effect from the effect 
of initial state and of the screening of the laser field and 
providing information on the factors involved. 
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